Quantum Phase Transitions of Hard-Core Bosons in Background Potentials 
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We study the zero temperature phase diagram of hard core bosons in two dimensions subjected 
to three types of background potentials: staggered, uniform, and random. In all three cases there 
is a quantum phase transition from a superfluid (at small potential) to a normal phase (at large 
potential), but with different universality classes. As expected, the staggered case belongs to the 
XY universality, while the uniform potential induces a mean field transition. The disorder driven 
transition is clearly different from both; in particular, we find z ~ 1.4, v ~ 1, and ji ~ 0.6. 
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The ground state of a quantum system can change 
abruptly at certain points when a parameter in the 
Hamiltonian is varied. Such "quantum phase transitions" 
(QPT) have attracted a great deal of attention recently 
[l| ; those that involve the interplay of disorder and inter- 
actions are among the most interesting and challenging. 

Superfluids and superconductors, for example, undergo 
a QPT to a normal state when the disorder in the 
medium is increased. Experimentally such superfluid- 
normal (superconductor-insulator) transitions occur in 
4 Hc films 0, arrays of Josephson junctions Q, and 
thin- film superconductors U The transition is typically 
induced by tuning the disorder strength, chemical po- 
tential, or magnetic field. Theoretically, it is generally 
accepted that this transition can be studied using the 
bosonic Hubbard model; techniques that have been used 
include mean field theory Jala, scaling and renormaliza- 
tion group argumen ts 13 , l a, Isl 1 1 Ct 1 1 ill - and Monte Carlo 
calculations [I^ll^ll4lllFlll6lll7l ll8 | . The transition is 
thought to be from the superfluid to a disordered "Bose 
glass" phase. 

Despite this extensive work, fundamental features of 
superconductor-insulator (SI) transitions due to disor- 
der, such as their university class (or classes), remain 
controversial. Of particular interest in a QPT is the dy- 
namical exponent z which relates the diverging correla- 
tion lengths in temporal and spatial directions: £t ~ £*. 
From scaling properties, under reasonable assumptions, 
it was argued that z = d, where d is the spatial dimension 
0- Those arguments, however, have been questioned re- 
cently 0, ED]- Early studies of the bosonic Hubbard 
model found conflicting results for the values of z and v 
[l2l, Il3| . More recently the "current-loop" version of this 
model has gained popularity 0, 0, 0, 0, 0] , espe- 
cially for investigating the effect of weak-disorder on the 
clean transition near integer filling [lH flfl flfl Il7j , and 
yields results generally consistent with those from scaling 
theory. 

In this paper we revisit the subject by providing a care- 
ful analysis of finite size scaling in and around the crit- 
ical region. In particular, compared to previous work, 
we study much bigger systems in the superfluid phase 



and make no assumption concerning z. The new method 
of analysis is first applied to two clean QPT in order to 
demonstrate its validity. We find that the strong-disorder 
SI transition that we study (hard core bosons at half fill- 
ing) has z « 1.4. This is strong evidence that disorder- 
induced SI transitions outside of the z = d universality 
class do exist, supporting 

Hard-core bosons in a background potential — We focus 
on a model of hard-core bosons in three different back- 
ground potentials. This is the boson Hubbard model 
in the limit of infinite repulsion and is equivalent to a 
current-loop model with restricted values of the current. 
The Hamiltonian is 
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where a] and dj are boson creation and annihilation op- 
erators for the site i = (i x , i y ) on a L x L bipartite square 
lattice and n, = aja^ is the number operator. The hard 
core nature of the bosons constrains the boson number at 
each site to be or 1. (Note that this model is equivalent 
to a spin-1/2 quantum XY model in a transverse mag- 
netic field.) We consider three types of background po- 
tentials (if. (1) a uniform potential yn = [i for all i, (2) a 
staggered potential given by /Lt« = o"i/x where cr, = 1 (or —1) 
when i belongs to the even (odd) sublattice, and (3) a dis- 
ordered potential given by [ii uniformly distributed be- 
tween —/it and fi. Note that in the case of a staggered 
potential, the system is exactly at half filling for all val- 
ues of /z; for the disordered potential, the system is at 
half filling on average but not for a given realization of 
the disorder (except in the thermodynamic limit). 

The above model has been studied with a uniform 
chemic al p otential [lflj . as well as with a staggered po- 
tential 2(;.i2JJ. However, the transition due to a random 
potential remains relatively unexplored by exact calcula- 
tions. The destruction of superfluidity by disorder was 
indeed studied in 0, but the universality of the 
phase transition was not addressed. A nice feature of 
our study is that we focus on three different universality 
classes in the context of a single simple model. 
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FIG. 1: The compressibility as a function of the magnitude 
of the three different background potentials. The uniform 
(circles) and staggered (squares) cases are very different from 
the average behavior in the disordered system (diamonds). 

We use the dircctcd-loop Monte Carlo algorithm in the 
discrete-time path-integral formulation [23j to compute 
quantities. In particular we divide the Euclidean time 
extent 1/T into M slices such that eM = 1/T. We fix 
e = 0.1 in this work. The partition function can then be 
written as a sum over configurations on an L x L x 4M 
lattice in the occupation basis. The extra four time slices 
arise when only nearest neighbor interactions are allowed 
between two consecutive time slices. 

The different behavior caused by the three background 
potentials is strikingly illustrated by looking at the com- 
pressibility as the magnitude of the potential is increased 
(see Fig.^). The compressibility is evaluated from 

i i 

The SI transition is very clear in the case of the ordered 
potentials: it occurs at p = 2 for the uniform field and 
near /z ~ 1 in the staggered case. At these points k 
becomes as the system develops a gap to excitations. 
In contrast, in the disordered case, no gap appears near 
the transition point, about 4.4, and so the compressibility 
remains non-zero. A compressible insulating phase is the 
signature of a Bose glass. 

Observables and scaling — We explore the physics of 
the model using three observables: (1) The winding num- 
ber susceptibility is 

X W = \(W 2 X +W*) (3) 

where W x and W y are the number of bosons crossing a 
surface perpendicular to the x and y directions respec- 
tively. In the thermodynamic limit this quantity is re- 
lated to the superfluid density by p s = T X w (2) The 
order-parameter susceptibility is 

Xp= m?M £ cfr'ni.'O ( 4 ) 

i,j,T,r' 




FIG. 2: Finite size scaling of \v an d X™ m the- superfluid 
phase for p close to p c in a staggered, uniform, and random 
potential (top to bottom). The solid lines show the fits to 
extract E 2 and p s . 

where C(i,T;j,r') — (ai(r)at(r')) is the boson correla- 
tion function between two space-time points, (i,r) and 
(j, r'). The superfluid condensate squared is defined by 
£ 2 = Txp/L 2 in the thermodynamic limit. (3) The spa- 
tial correlation length £l is defined using the relation 

G(i x ,j x ) = C ^ ^ 3, r') oc cosh ( |L/2 = ~ jx) 1 ) 

■i t' 

(5) 

for large spatial box size L and separations \i x — j x \. We 
refer to the correlation length at L — oo as simply £ and 
compute it through a large L extrapolation of £l . 

In order to access T = properties, we study the be- 
havior of Xwi Xpi an d £ as a function of L, T, and fi. In 
the disordered case, we average over at least 50 disorder 
realizations. (For smaller lattices, we use several 100 re- 
alizations.) To check that this number is sufficient, for a 
few values of \i close to the transition we used up to 5000 
realizations. In all these cases we found that the only 
role of increasing the number of disorder realizations in 
computing the average is to reduce the error bars, which 
shrink as expected, i.e., as 1/yN for N > 40 realizations. 

For a second order phase transition, in the critical re- 
gion (5 = (/x — jtt c ) ~ 0, T — * 0, L —>■ oo) scaling theory 
predicts |l| 

Xw = fiSL 1 ^, TL Z ), Xp = £ 2 -". 9 («^, TL Z ). (6) 

When 8 < 0, there is a critical temperature for su- 
perfluidity such that for T < T c (8) and as L — > oo, 
Xw ~ Ps/T and Xp ~ S 2 L 2 jT . Consistency then con- 
strains the scaling functions: ]{8L 1 I V \TL Z ) - \8\ zv jT 
&ndg(8L 1 / v ,TL z ) - \8\ 2f3 /T where 2/3 = v(z + rj). Thus, 
we expect p s = A|J| zly and S 2 = B^l 2 ' 5 . In order to ex- 
tract p s and S 2 , we can set 1/T = L and study the be- 
havior of Xw and Xp as a function of L. Note that in the 
superfluid phase we do not have to set 1/T ~ L z , making 
possible calculations up to L = 96. We fit our data to the 
form Xw — Lp s +a + b/L and Xp = L 3 T, 2 + cL? + dL for 
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FIG. 3: Uniform potential: p s , E 2 , and £ _1 as a function of 
the potential strength. The QPT is at p C:U = 2.0. The inset 
shows the scaling of the correlation length with temperature 
at p c ,u- The solid lines are a joint fit to the expected mean 
field scaling. 

large L. Our data fits well to this expectation, allowing 
p s and E 2 to be reliably extracted. Figure [21 shows this 
finite size scaling behavior for each of the three types of 
potentials at a value of p such that \p — p c \/ 1 p c < 5%. 

When (5 = scaling theory predicts that Xw is a univer- 
sal function of a single variable TL Z with \w = fo/(TL z ) 
for TL Z < 1 and X w = fi exp(-aT 1 / z i) for TL Z > 1. 
We see that at 8 = the correlation length £ in the 
L — > oo limit is controlled by the temperature and di- 
verges according to = CT 1 ^ as T — > 0. We find z 
by first extracting £ by using calculations up to L = 96 
and then fitting £ to the asymptotic form above. When 
S > 0, £ is expected to be finite even at T = 0. In that 
regime, the relation = Do~ v defines the exponent v. 

Uniform background potential — The QPT driven by 
a uniform chemical potential in a system of hard core 
bosons is expected to be described by mean field the- 
ory in two or more spatial dimensions 0. Since two 
spatial dimensions is the upper critical dimension, one 
can have large logarithmic corrections to the usual scal- 
ing forms. We assume p s = A\6\ z "(l + A'/logS), E 2 = 
B\5\ 2I3 {1 + B'/logS), = CT Y I Z (\ - C'/logT) and 
£" 1 = D5 V (1 + D'/\og5), where i/ = 1/2, = 1/2 and 
2 = 2. Further, it is known (lij that p c — 2. Our data, 
shown in Fig. |3| fits these expectations very well in the 
range 1.85 < p < 2.1, L > 6 (x 2 /DOF = 0.5) 0. 

Staggered background potential — A system of hard core 
bosons in a staggered chemical potential has a remain- 
ing Zi symmetry (translation by one unit followed by a 
particle- hole transformation). Due to this symmetry, the 
QPT is expected to be in the three-dimensional (3d) XY 
universality class with v — 0.672, (3 = 0.348, and z = 1. 
With these fixed exponent values, a combined fit of all 
our data (exce pt f or L < 4), shown in Fig. works 
extremely well 25j. Note that the data to demonstrate 
£ = CT X I Z (inset of Fig. was obtained slightly off p c ; 
we find that this does not influence the value of z. 




FIG. 4: Staggered potential: p s , E 2 , and £ 1 as a function 
of the potential strength. The QPT is at p c>s = 0.9919(4). 
The inset shows the scaling of the correlation length with 
temperature at p — 0.991. The solid lines are a joint fit to 
the expected 3dXY scaling. 

Disordered potential — Assuming that the universality 
class is unknown for a disordered background potential, 
we fit our data for p s , E 2 , and £, shown in Fig. \5\ to the 
expected scaling forms. Note that £(/Lt) is missing: £(T) 
appears not to saturate at small T for p > /i c , and thus 
we were unable to extract £(T = 0) reliably. 

Interestingly, we can fit the data for a large range, 
p > 3.3, to the scaling forms: a combined fit yields 
p c = 4.42(2), v = 1.10(4), (3 = 0.61(2), and z = 1.40(2) 
with x 2 /DOF = 0.85 (see for prefactors). Since 
2(3 = v(z + rj), we estimate rj ~ —0.22(6). Note that 
z is extracted by fitting ^ = CT 1/z at p = 4.42. Al- 
though p c may be slightly different, we believe the value 
of z should not be affected much, based on our experi- 
ence with the staggered case. Further support for this 
value of z is shown in Fig. [S] at p c > one expects XwTL z 
to be a function of one variable, TL Z , and we see a good 




FIG. 5: Disordered potential: p s and E 2 as a function of the 
potential strength. The QPT is at p c , d = 4.42(2). The inset 
shows the scaling of the correlation length with temperature 
at p = 4.42. The solid lines are a joint fit to critical behavior; 
the dashed line is a fit to p s = A(p Ct d — p) vz assuming v = 1 
and z — 2. 
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FIG. 6: Scaling function for Xw sA, fj, = fj, c choosing z = 1.4 
(main figure) or z = 2 (left inset). The right inset shows 
the decay of T\ m when T is scaled as 2/L 2 ; the solid line is 
0.93/L 1 ' 4 . 

collapse of the data when we use fi = 4.42 and z — 1.4. 

Clearly, the universality class of the phase transitions 
with staggered, uniform, and random potentials are strik- 
ingly different. The large size of the critical region for a 
disordered potential is a bit surprising. Note, however, 
that our results are consistent with the expectation that 
v > 2/d ■ On the other hand, we do not see z = 2 as 
predicted by 0. Indeed, for 2 = 2, \wTL z is no longer a 
function of TL Z (inset to Fig.EJ). Earlier work for strong 
disorder using a current loop model found z ~ 2 and 
v ~ 1 18], implying p s — A\S\ 2 . We can fit our data 
to this form (dashed line in Fig. [SJ) if we narrow the 
fitting region to fi > 4. Such a fit yields /i c = 4.49(2); 
This would imply that if we scale according to T = 2/L 2 
then at [i = 4.42 we must find T\w goes to a constant 
for large L. However, we find that T\w ~ 1/L 1A (in- 
set to Fig. 01 for lattices as large as L = 32 (for which 
1/T = 512), suggesting that fi c is close to 4.42. In fact, 
even at fi = 4.6, T Xw ~ 1/L 1M for T = 2/L 2 scaling. If 
we make the assumption that \i c < 4.6, then we get the 
bound z < 1.6. 

Three possible explanations come to mind for the dif- 
ference between our results and other recent results: 
(1) A fixed choice of z introduced a bias in the previous 
analysis. In particular, a combination of small lattices 
and a shifted critical point may have conspired to give 
results consistent with z — 2. (2) The disorder averages 
which we have performed may be insufficient: there may 
exist statistically significant rare events which we have 
not encountered. (3) The real critical region may be en- 
countered only for /i > 4.35 for which we would need 
lattices much larger than is possible to date. 
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